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FOREWORD 


The work reported herein was performed by the Grumman Aerospace Corporation 
under the NASA/Langley Master Agreement and Contract No. NAS 1-10635 for the 
Development and Implementation of Space Shuttle Structural Dynamics Modeling 
Technology. The Work Statements of Task Order No. 17, "Development of An 
Analytical Program to Analyze Reusable Surface Insulation for Shuttle," and Task 
Order No. 19, "Development of a Thermal Stress Analysis Program for Reusable 
Surface Insulation for Shuttle", authorized and specified the tasks to be per- 
formed in this study. The period of performance for this effort was for 15 
months starting in June of 1973* 

The overall supervision of programs under the Master Agreement is provided 
by Mr. E. F. Baird, Master Agreement Program Manager. The Project Manager for 
Task Order Nos. 17 and 19 was Dr. I. U. Ojalvo. Many individuals at Grumman con- 
tributed to the work reported herein. However, the authors wish to specifically 
acknowledge the Applied Mechanics Group of the Grumman Research Department for 
supplying some of the basic finite element subroutines required for this work. 
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I. INTRODUCTION 


This manual presents programming information for the RESIST (REusable 
Surface Insulation STresses) program for the dynamic and thermal stress analysis 
of the Space Shuttle surface insulation. A user's manual which describes the 
necessary input is contained as Appendix B in both References 1 and 2. The engi- 
neering solution procedure programmed is also described in these reports. 

The present manual contains the overall flow chart of the program, overlay 
chart, data set allocation and subprogram calling sequence, along with a brief 
description of the individual subprograms and typical subprogram output. Two 
versions of the program exist, one of which is written completely in FORTRAN. 

This version is compatible with both the IBM 370/l68 and CDC 6600 machines. A 
more efficient second version exists which runs only on the IBM computer since it 
contains several assembly language subroutines. These routines are: 

• DINIT with entry points: DREAD 

PREAD 

DWRITE 

PRITE 

REND 

WEND 

DFIND 

DCLOSE 

• SECOND with entry point ICHRON 
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II 


FLOW CHART 


The overall logic for the RESIST computer program is displayed schematically 
through use of 12 modules and 6 logical and arithmetic "if" (diamond- shaped) boxes 
in Figure 2-1. Input data for Module I is detailed in References 1 & 2. Finite 
element topological information for the RSI tiles and primary structure is auto- 
matically generated in Module II. The element stiffness matrices for the primary 
structure and a typical RSI tile are generated and assembled in Module III, after 
which a test is performed to establish if the statics (ST) or modal vibration (VB) 
option is to be exercised. * 

When the vibration option is operational, the index (designated NH in the 
program) is set to zero and the primary structure mass matrix and two tile mass 
matrices are generated. The primary structure mass matrix is then combined with 
the tile mass matrix, which only approximates, tile inertia effects. The approxi- 
mation used is associated with the neglect of tile stiffness and a kinematic 
assumption, which is tantamount to assuming that the tile deformations are devoid 
of shear effects for the purpose of computing their modal kinetic energy. Approx- 
imate mode shapes and frequencies for the primary structure are then computed in 
Module V (designated ALARM and detailed in Reference l). 

The iteration index is then stepped up by 1 and the primary structure deflec- 
tions are imposed at the tile/primary-structure interfaces, and the tile deflections 
are obtained in Module VI. This procedure is performed for all tiles, or only for 
those designated by an input option. If all tiles are treated in Module VI, a 
Rayleigh Quotient is computed next (Module VII ) to update the system frequency. 

A tile- loading matrix is then formed (Module VIII ) and applied to obtain new 
primary structure deflections in Module IX. 

Primary structure deflection differences from the previous iterate’ s results 
are checked for convergence against an input tolerance. If convergence has been 
achieved, the program computes tile and plate stresses when these are requested 
by the user. If convergence is not achieved, a test is performed to determine 
whether the maximum number of iterations have been performed. The program recycles 
through Modules VI through IX until convergence is achieved or the maximum number 
of iterations have been performed. 


When the statics option is selected, the index NH is set to 1 after Module III 
The primary structure loading is then computed. in Modules X and VIII. The proced- 
ural logic then cycles through Modules IX, VI and VIII until convergence or the 
maximum number of iterations is obtained. Regardless of which type problem is 
solved, ST or VB, all program modules with the exception of IV, V, VII and X are 
shared. 
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Figure 2-1 Flow Chart for RSI Stress Analysis Program ''RESIST' 
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III. OVERLAY CHART 

The following is the overlay chart for the RESIST program as used on an 
IBM 370/168. This chart contains the root segment plus five levels of overlay 
(denoted by A through E), and depicts the subprograms and labeled common blocks 
(enclosed within slashes, e.g. /KTAB,/) of each segment, along with the starting 
address and length (in bytes) of each segment. The numbers in the upper right 
hand corner of each box denote the segment number. This overlay structure also 
applies for the CDC 6600 except that the segment lengths are different. 
























IV. DATA SET ALLOCATIONS 

The following chart identifies the unit and file number assignments for 
each data set. Symbols used in the chart are: 

BC Boundary Condition 

DOF Degrees of Freedom 

K Stiffness Matrix 

M Mass Matrix 

P Load Vector 

PS Primary Structure 

6 Deflections 

u) Natural Frequency 
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FILE 


UNIT 

UNIT 


_2_ 

_3_ 

NO. 

NAME 




1 

JDSRN 1 

PS boundary conditions (No. nodes x 1 1 ) 



2 

ISCR 

Normalized PS mode shapes (No. modes x No. DOF) 

Partition of tile K— o> 2 M into K AA 
/ No. tile DOF — No. bound DOF x \ 
\No. tile DOF — No. bound DOF J 


3 

JSCR 

PS load vector with applied BC's (No. DOF x 1) 



4 

MDSRN 

PS loads with PZ = 0 before BC's are applied 

PS total stiffness matrix with zeroes in upper 




(No. nodes x 6) 

triangle (No. DOF x No. DOF) 


6 

ITAPEW 

Output to printer 



7 

NPIRI 

Nodes, coordinates & matl properties for each tile 
member (No. members x 81) 



8 

NDSRN 

1 PS element mass matrices (No. elements - No. rows) 





2 PS total diagonal mass matrix (No. DOF x 1) 



9 

LDSRN 

1 PS element stiffness matrices (No. elements = No. rows) 

Scratch 




2 PS total mass matrix (No. DOF x No. DOF) 



10 

IDSRN 

1 PS geometry (No. nodes x 5) 

2 Decomposed lower triangle of PS 





K-Matrix (No. DOF x No. DOF) 

PS deflections (No. DOF x 1) 


11 

IUNIT 

Partition of tile K— oj 2 M j nto 
Kab/No. tile DOF — No. bound DOF x \ 

Partition of tile K— u> 2 M 




yNo. bound DOF J 

into K bb 

/ No. bound DOF x \ 
\^No. bound DOF J 


12 

JUNIT 

1 Tile element stiffness matrices 
(No. elements = No. rows) 




NPIRF 

2 Scratch 



13 

KUNIT 

Full tile K-matrix / No. tile DOF x \ 

Partition of tile load vector into P A 

Partition of tile load vector 



\No. tile DOF ) 

^No. tile DOF — No. bound DOF x 1^ 

into P B 

(no. bound DOF x 1 ^ 

14 

LUNIT 

1 For tile: K aa [Pa~K ab 8 b ] 





/ No. tile DOF — No. bound DOF x \ 
\No. tile DOF — No. bound DOF J 




NPIRK 

2 PS deflections associated with individual tile (No. 

Boundary loads generated in tile 




tile bound DOF x 1 ) 

analysis (No. bound DOF x 1) 


15 

ITAPER 

1 Input data set in card image format (20A4) 




NPIRA 

2 Tile load vector (No. tile DOF x 1) 



16 

NTAPEI 

Input stream 



17 

NPIRB 

1 PS members (No. elements x 100) 





2 Tile (K -0) 2 M) 





^No. tile DOF x No. tile DOF^) 



18 

NPIRC 

Scratch 



19 

NPIRD 

Tile displacements - boundary displacements 





^No. tile DOF - No. bound DOF x 1 ) 



20 ‘ 

NPIRE 

Tile strain matrix (No. elements = No. rows) 




NPIRG 


Tile total mass matrix 

Tile total stiffness matrix 


(1 x No. tile DOF) 


with zeroes in upper triangle 
(No. tile DOF x No. tile DOF 


Figure 4-1 Data Set Allocations 
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V. SUBPROGRAM CALLING SEQUENCE AND COMMON BLOCKS UTILIZED 

The following table lists the subprograms that are called by each parent 
routine. The common blocks used in each routine are also tabularized. The sub- 
programs listed under the "Calls" heading appear in the order in which they are 
called. The order of the subprograms listed in the "Subprogram" column corre- 
sponds to the order in which they appear in the FORTRAN listing. An alphabetical 
index which cross references this order-of-appearance number is provided in 
Section VI. 
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SUBPROGRAM 

CALLED FROM 


CALLS 

COMMON 

MAIN 

- - 


TIME A' " 

CTAPES 




BIGBRD 

CTMH 




LDB 

CLIST 




EOFOl 

C ID I V 




SAMAIN 

CBYTES 

CORE 





UNCUTT 

S AMAIN 

main 

“ 

SETUP 

FILE 




CARDIN 

MATRIX 




MESH 

CTAPES 

*■ - ' 



ELAS 
"TIME A' 
PROCES 
INPUT 
SBMAIN 

RWO 

SETUP 

samain 



DSRN 

UNITS 

meminf 

FILE 

MATRIX 

CARDIN 

SAMAIN 


mprop 

CTAPES 




GEOBC 

DIMEN 




MEMBIN 

MATRIX 




LOADIN 

UNCUTT 

INFO 

AFLEX 

PR0BS2 

RWO 

GEOBC 

CARDIN 


BRICK 

FILE 




ijndcut 

CTAPES 




putlab 

SIZE 




putrow 

KTAB 




DCLOSE 

MINBOF 




BOUND 

DIMEN 

MATRIX 

PROBSZ 

BRICKT 

UNCUTT 

BRICK 

GEOBC 



MINBOF 

BRICKT 

DIMEN 

PROBSZ 

CTAPES 

UNOCUT 

GEOBC 



MINBOF 

UNCUTT 

DIMEN 

PROBSZ 

CTAPES 
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NO. 

SUBPROGRAM 

CALLED FROM 

CALLS 

COMMON 

8 

BOUND 

GEOBC 

.... 

CTAPES 


SIZE 

DIMEN 

PROBSZ 

KTAB 

AFLEX 


9 

MEMBIN 

CARDIN 

PUTLAB 

FILE 




MEMGEN 

MATRIX 




DCLOSE 

PROBSZ 

SIZE 

DIMEN 

10 

MEMGEN 

MEMBIN 

PUTROW 

MATRIX 




SORT 

PROBSZ 





MINBOF 

INFO 

DIMEN 

11 

LOADIN 

CARDIN 

PUTLAB 

FILE 




PUTROW 

CTAPES 




DCLOSE 

MINBOF 

SIZE 

PROBSZ 

CBYTES 

DIMEN 

INFO 

RWO 

MATRIX 

KTAB 

12 

MESH 

SAMAIN 

INTEMP 

DIMEN 




FTEMP 

INFO 




PUTLAB 

UNCUTT 




MSHWRT 

UNITS 




DCLOSE 

CTAPES 

MEMINF 

PANDR 

CBYTES 

LENGTH 

13 

MSHV.'RT 

MESH 

TINIER 

CTAPES 




matl 

MEMINF 




PUTROW 

PANDR 

TABIL 

INFO 

14 . 

MPROP 

CARDIN 


INFO 

TABIL 

DIMEN 

RWO 

CTAPES 
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""NO. 

SUBPROGRAM 

CALLED FROM 

CALLS 

COMMON' 

--re 

— MATL 




* MSHWRT" 



16 

TIMER 

MSHWRT 

LAGRAN 

TABIL 

CTAPES 

17 

LAGRAN 

tinter 

— 

— 

18 

INTEMP 

MESH 


TTAB 

CTAPES 

CELAS 

19 

FTEMP 

MESH 

TEMPER 

TTAB 





CTAPES 

20 

TEMPER 

FTEMP 

TLGRN 

TTAB 

21 

TLGRN 

TEMPER 

DLGRN 

TTAB 




LGRANG 


22 

DLGRN 

TLGRN 

LGRANG 

— 

23 

LGRANG 

TLGRN 

DLGRN 1 

— 

— 

- 24 

COPY 

SBMAIN 

GETDIM 




HEXTIL 

PUTLAB 

GETROW 

PUTROW 

DCLOSE 


23 

ELAS 

SAVA IN 

SLOAD 

MEMINF 


- 


fsakg 

UNITS 




PARTIN 

PRINT 


OFACT CELAS 

LENGTH 

INFO 

PANDR 

PROBSZ 

FILE 

RWO 


26 

FMSTS 

HEXEL 

JACOB 

-• — - 

27 

FSAKG 

ELAS 

HEXEL 

MEMINF 




PUTLAB 

UNITS 




PUTROW 

PRINT 




DCLOSE 

NSTAK 

RESTOR 

file 


28 GAUSS 


HEXEL 


C TAPES 


5 - 1 + 


NO. 

SUBPROGRAM 

CALLED FROM 

CALLS 

COMMON 

29 

HEXEL 

F5AKG 

PUTLAB 

INFO 




GETDIM 

MEMINF 




GETROW 

UNITS 




GAUSS 

CELAS 




SHAPE 

JACOB 

ODST 

ODLD 

WRAKG 

ISTAK 

FMSTS 

PRITE 

WEND 

DOLOSE 

PANDR 

30 

ISTAK 

HEXEL 

MSTAK 

MEMINF 




PRITE 

PRINT 




. WEND 

UNITS 

CELAS 

CTAPES 


31 

JACOB 

HEXEL 

FMSTS 



32 

MSTAK 

ISTAK 

— 

PRINT 

33 

NSTAK 

FSAKG 

PUTLAB 

GETDIM 

PREAD 

REND 

PACK 

PUTROW 

DOLOSE 

UNITS 

PRINT 

3 A 

PARTIN 

ELAS 

HEXTIL 

GETDIM 

PUTLAB 

GETROW 

PUTROW 

DOLOSE 


35 

ODLD 

HEXEL 

— 


36 

QDST 

HEXEL 

— 

— 

37 

SHAPE 

HEXEL 

— 


38 

SLOAD 

ELAS 

— — 

MEMINF 

CELAS 

PRINT 





CTAPES 

39 

WRAKG 

HEXEL 

— — 

PRINT 


CT APES 
MEMINF 
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(TO. SUBPROGRAM CALC ED FROM” 1 CALCS'"' 


COMMON 


4U~ 

— P RuCES 

SANA IN 

GETDIM 

GETROW 

TITLES 

JTAB 

FILE 

MATRIX 

trobsz - 

CTAPES 

CLIST 





41 

INPUT 

SAMAIN 

GETDIM 

PUTLAB 

“* GETROW 

TITLES 

EL2 

JTAB 

FILE 

MINBUF 

MATRIX 

PROBSZ 




“EU5" . ' 

DIMEN 




EL 16 

INFO 




PLB 

CTAPES 




DC LOSE 

' CLIST 

42 

SPLITS 

EL2 

FLXIBL 

JTAB 



“EL5~ 

' A RACE 

'MTNBUF" 



EL16 

EDGMAS 

MATRIX 




PRINT3 

PROBSZ 




PRITE 

INFO 




WEND 

CTAPES 





CBYTES 

43 

EL2 

INPUT 

DSINF 

MINBUF 




DCOSF 

PROBSZ 




DSQRTF 

CTAPES 




SPLITS 

INFO 

44 - 

• ” EL 5 - 

INPUT 

DSINF 

MINBUF 


- 


DABSF 

PROBSZ 




DCOSF 

CTAPES 




DSQRTF 

INFO 




SPLITS 


' 45 - 

ELI 5 - * • 

ELI 6 

DSQRTF 

PROBSZ 





MINBUF 





ARRAYS 





STIFF 





CTAPES 

46 

-ELI6 

INPUT 

DSQRTF 

PROBSZ 




DABSF 

MINBUF 




DCOSF 

ARRAYS 




DSINF 

STIFF 




EL 1 5 

CTAPES 




SREVN2 

INFO 




SPLITS 


47 

SREVN2 

EL16 

DABSF 

— 

48 

FLXIBL 

SPLITS 


AFLEX 





DIMEN 


\ 
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SUBPROGRAM 


CALLED FROM 


CALLS 


COMMON 


NO. 


49 

EDGMAS 

SPLITS 

““ 

AFLEX 

50 

ARACE 

SPLITS 

— 

JTAB 

51 

PRINT3 

SPLITS 

PLB 

CT APES 




TITLES 

CLIST 

52 

SBMAIN 

SAMAIN 

ASTACK 

DSRN 




TIMEA 

FILE 




QFACT 

MATRIX 




MSTACK 

UNITS 




MPR I NT 

PROBS2 




BPMASS 

UNCUTT 




ALARM 

INFO 




RESTOR 

R WO 




GETDIM 

C TAPES 




PUTLAB 

GETROW 

CORE 


- • 


PUTROW 

DCLOSE 

LOCDEF 

MULT 

READMT 

ICMULT 

QFACT 

KMOSQ 

COPY 

QFSOL 

REVERS 

QBSOL 

ORTHOG 

MSOUT 

COMPAR 

PSTRES 

SORT 


53 

ASTACK 

SBMAIN 

PUTLAB 

PROBS2 




GETDIM 

DIMEN 




PREAD 

INFO 




REND 

C TAPES 




PACK 

PUTROW 

DCLOSE 

CBYTES 

54 

MSTACK 

SBMAIN 

PUTLAB 

PROBSZ 




GETD I w 

C TAPES 




PREAD 

REND 

PACK 

PUTROW 

DCLOSE 

CBYTES 
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TCOT SUBPROGRAM - ' CALTTFD 'FROM CALLS ' COMMON" 


55 

- *\ 

"BPMA55 ' 

SBMA I N 


GETDIM . ' 

PREAO 

REND 

. PROBSZ 
CBYTES 





potlaB 

PUTROW 

DCLOSE 


56 

ALARM 

sbmain 


TITLES 

REVERS 

TRAN" 

PUTLAB 

SCAPRO 

PROBSZ 

INFO 

FTLT 

CTAPES 

CLIST 





5QR T 
PUTROW 
DCLOSE 
QFSOL 






GETOIM 

GETROW 

MULT *“* ' " 






ABS 
AMAX1 
AM INI 






STURM 

QSVEC2 

MSOUT 

_ 

57 

SCAPRO 

ALARM 


ROM 


58 

STURM 

ALARM 


ABS 

— — mm— 





PREP 


59 

PREP 

STURM 


— 

— 

60 

QSVEC2 

ALARM 


ABS 

PATLVN 

— - — • — 




amaxi 

SORT 

• - — 

61 

DOTPRO 

QSVEC2 


RDM 

ANOOR 

DOTPRO 

PATLVN 

62 

ANOOR 

QSVEC2 


— 


-• 63 

ROM 

ALARM 


__ _ _ 


6r4 

RESTOR 

0SVEC2 

FSAKG 


GETDtM ~ 

CTAPES " 

•• 


SBMAIN 


GETROW 

PUTLAB 

PUTROW 

DCLOSE 
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SUBPROGRAM 


CALLED FROM 


CALLS 


COMMON 


NO. 


65 MULT 


SBMAIN 

ALARM 

HEXTIL 


GETDIM 

GETROW 

PUTLAB 

ENMMPY 

PUTROW 

DCLOSE 


CTAPES 


66 ENMMPY MULT 

67 READMT SBMAIN GETDIM 

GETROW 

DCLOSE 


68 ICMULT 


SBMAIN 

HEXTIL 


69 LOWTRI 


SBMAIN 

HEXTIL 


GETDIM 

PUTLAB 

GETROW 

PUTROW 

DCLOSE 


70 

QFACT 

ELAS 

GETDIM 

CTAPES 



SBMAIN 

PUTLAB 

CIDI V 



PODSYM 

OCHOL 

INFO 




DCLOSE 


71 

OCHOL 

OFACT 

TITLES 

PROBSZ 




GETROW 

CTAPES 




UNPACK 

R WO 




KPRINT 

CLIST 




HOTDOT 

INFO 




PUTROW 





ALOGIO 





SORT 





ARS 





IABS 


72 

KPRINT 

OCHOL 


CTAPES 

73 

OPASS 

GFSOL 

GETROW 

— 

74 

OF SOL 

SBMAIN 

GETDIM 

FILE 



ALARM 

PUTLAB 

MATRIX 



PODSYM 

OPASS 

CTAPES 




DCLOSE 

C I D I V 


QFOR 
IABS 
MAXO 
MI NO 
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~Ncr. 

— 73 ~ 


76 


,77 

78 

79 

80 

81 


SUBPROGRAM CALLED FROM CALLS COMMON 


QFOR - - 

OFSOLT 

GEJRm - 

UNPACK 

HOTDOT 

PUTROW 

MAXO 

MINO 


QBSOL 

SBMAIN 

GETDIM 

CTAPES 


PODSYM 

PUTLAB 

C J D I V 



QBAC" 




DCLOSE 




I ABS 


r _ ~ 


MAX'D 


QBAC 

QBSOL 

GETROW 



UNPACK 

HOTDOT 

PUTROW 


REVERS 

SBMAIN 

GETDIM 



■ ALARM 

PUTLAB 



PODSYM 

GETROW 

PUTROW 




DCLOSE 


TRAN • 

ALARM 

GETDIM 

GETROW 

PUTLAB 

PUTROW 




DCLOSE 


ORTHOG 

SBMAIN 

GETDIM 

GETROW 

SORT 

PUTLAB 

PUTROW 




DCLOSE 


LOCDEF 

SBMAIN 

GETDIM 

DIMEN 



PUTLAB 

FILE 



GETROW 

CTAPES 



PUTROW 

CLIST 



DCLOSE 

MATRIX 



DWRITE 

UNITS 



HEXTIL 

RWO 



DREAD 

INFO 


subprogram 


CALLED FROM 


CALLS 


COMMON 


NO. 


a 2 

hextil 

LOCDEF 

KMOSQ 

MEM INF 




PARTIN 

UNITS 




COPY 

PRINT 




MULT 

PANDR 




LOW TR I 

INFO 




SUB 

RWO 

■- 



PODSYM 

CTAPES 




MATB 

AOO 

GETDIM 

GETROW 

DCLOSE 

FILE 




WRDSP 
ICMULT 
SR A I N2 
SRAIN 
SRAIN3 


83 

KMOSO 

SBMA IN 

GETDIM 




HEXTIL 

GETROW 

putlab 

PUTROW 

DCLOSE 


84 

OSRN 

SRAIN 

— 

— 

8$ 

srain 

HEX TIL 

GETDIM 

UNITS 




GETROW 

CTAPES 




DREAD 

CELAS 




PREAD 

REND 

DSRAN 

DCLOSE 

RWO 

06 

SRAIN2 

HEXTIL 

G E T D I M 

UNITS 




GETROW 

CTAPES 

CELAS 

RWO 

87 

SRAIN3 

HEXTIL 

GETDIM 

UNITS 




GETROW 

CTAPES 




DCLOSE 

CELAS 

RWO 

88 

WRDSP 

HEXTIL 


MEM INF 
RWO 

CTAPES 

89 

PODSYM 

HEXTIL 

QFACT 

REVERS 

QFSOL 

OBSOL 
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"“NO* SUBPROGRAM CALCED FROM 

KExrrt — 


CALLS'” 


"COMMON 


"GETVTVT 

PUTLAB 

GETROW 

PUTROW 

DCLOSE 


95 


96 


97 


98 


SUB 

HEXTIL 

. GET-DI-M 

PUTLAB 

GETROW 

PUTROW 

DCLOSE 



MATB“ 

HEXTIL ' 

GETD I M 
GETROW 
PUTLAB 
ENTMPY 
PUTROW 
DCLOSE 



entmpy 

MATB 

— 




PSTRES 

SBMAIN 

GETDIM 

GETROW 

DCLOSE 


DIMEN 
FILE 
C T APES 





RWO ~ "" 

PROBS2 

MATRIX 

MSOUT 

SBMAIN 

ALARM 

GETDIM 

GETROW 


MATRIX 

DIMEN 



DCLOSE 


CTAPES 

RWO 

FILE 

COMPAR 

SBMAIN 

GETDIM 

GETROW 


CTAPES 



DCLOSE 

ABS 



MPRINT 

SBMAIN 

TITLES 

GETDIM 

GETROW 


CTAPES 

CLIST 



UNPACK 

DCLOSE 



HOTOOT 

QCHOL 

QFOR 

OBAC 




— 
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NO •" ' SUBPROGRAM CALLED FROM 


CALLS 


COMMON 


PUTLAB ' 

PARTIN' 

DFTND' 


“PUT GET 


MULT 

DWRITE 


FILE 


COPY 



PROBS2 


GEOBC 



CTAPES 


LOAD IN 

SBMAIN 

QFSOL 

REVERS 

KMOSO 

ALARM 

FSAKG 

ASTACK 

NSTAK 

ADD 

MATB 

RESTOR 

QFACT 

INPUT 

MEMBIN 

MESH 

LOWTRI 

QBSOL 

TRAN 

ORTHOG 

BPMASS 

MSTACK 

LOCDEF 

HEXEL 

SUB 
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NO." SUBPROGRAM CALLED FROM 


CALL'S 


COMMON 


TOO - — PUTROW " PARTIN' DWRI7E 

MULT PACK 

COPY UNPACK 

mEMGEN DCLOSE 

MSHWRT JABS 


LOWTRI 

OB AC 

TRAN 

ORTHOG 

BPMASS 

MSTACK 

LOCOEF 

ADD " 

MATB 

RESTOR 

OCHOL 

GEOBC 

LOADIN 

SBMAIN 

OFOR 

REVERS 

K M OSO 

ALARM 

FSAKG 

ASTACK 

NSTAK 

SUB 
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NO* 

SUBPROGRAM 

CALLED FROM 

CALLS 

COMMON 

101 

GETDIM 

PARTIN 

OFIND 

PUTGET 



MULT 

DREAD 

FILE 



COPY 


PR0BS2 



INPUT 


CT APES 


LOWTRI 

OBSOL 

TRAN 


READMT 

PSTRES 

ALARM 

MSTACK 

LOCDEF 

HEXEL 

hextiL 

ADD 

SUB 

SRAIN2 

RESTOR 

QFACT 

PROCES 

SBMA I N 

QFSOL 

REVERS 

KMOSO 

ORTHOG 

MPRINT 

bpmass 

ASTACK 

NSTAK 

MSOUT 

COMPAR 

SRAIN 

MATB 

SRAIN3 


5-15 



~NOi SUBPROGRAM CALLED FROM 


CALLS 


COMMON 


T02' GETROW"^ "" PARTIN DREAD 

MULT UNPACK 

COPY PACK 

INPUT IABS 


LOWTRI 

QFOR 

REVERS "* - ’ 

KMOSO f 

ORTHOG 

— M PRINT 

LOCDEF 

MSOUT 

---- COMPAR 

SRAIN 

MATB 

SRAIN3 

RESTOR 

OCHOL 

" ; PROCES 

SBMAIN 
OP ASS 
OBAC 
TRAN 
READMT 

PSTRES 

ALARM 

HEXEL 

HEXTIL 

ADD 

SUR 

5RAIN2 

103 PRITE SPLITS PUT 

HEXEL 

ISTAK 

104 DWRITE PUTLAB PUT 

PUTROW 

LOCDEF 


105 PUT 


PRITE 

DWRITE 


106 PREAD MSTAK GET 

NSTAK 

BPMASS 

ASTACK 

SRAIN 


107 DREAD GETDIM GET 

GETROW 

LOCDEF 

SRAIN 
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NO# 

SUBPROGRAM 

CALLED FROM 

CALLS 

COMMON 

108 

GET 

PREAD 

DREAD 

— 


109 

WEND 

SPLITS 

HEXEL 

ISTAK 

.... 

.... 

110 

REND 

MSTAK 

NSTAK 

BPMASS 

ASTACK 

SRAIN 



111 

DFIND 

PUTLAB 

GETDIM 



— 

112 

DOLOSE 

PARTIN 

MULT 

— 

— 


PUTROW 

INPUT 

MEMBIN 

MESH 

LOWTRI 

Q3S0L 

TRAN 

READMT 

PSTRES 

ALARM 

FSAKG 

ASTACK 

NSTAK 

MS OUT 

COmPAR 

SPAIN 

MATR 

RESTOR 

QFACT 

COPY 

GECBC 

LOAD IN 

SBMAIN 

OFSOL 

REVER5 

K^OSQ 

ORTHOG 

MPPINT 

PPMASS 

MSTACK 

LOCDEF 

HEXEL 

HEXTIL 

ADD 

SUB 

SRAIN3 
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NO. 
' "113 

I1A 


115 

116 

1 17 
118 

119 

120 

121 

122 

123 

124 


SUBPROGRAM 


CALLED FROM 


CALLS 


common 


PACK 

PUTROW 

GETROW 

ASTACK 

MSTACK 

NSTAK 



UNPACK 

PUTROW 

GETROW 

QCHOl 

MINO 

— 


- — QFOR — 

CRAC 

MPRINT 



LDB 

MAIN 

EOFOl 

PLB 

— 

PLB 

TITLES 

INPUT 

LDR 

TIMEA 

PRINT3 


.... 

EOFOl 

MAIN 
LDR . 

— 

— 

TITLES 

TIMEA 

INPUT 

MPRINT 

QCHOL 

PROCES 

PRINT3 

ALARM 

PLB 

ctmh 

CLIST 

CTAPES 

PIGPRD 

main 

— 

ctapes 

timea 

main 
s amain 
SR VAIN 

TITLES 

PLB 

C TAPES 

CLIST 

PR08SZ 

DARSF 

EL5 
EL 16 
SREVN? 

ARS 

— 

DCOSF 

EL2 

ELS 

EL16 

COS 

— 

DSINF 

. EL? 
EL5 
EL 16 

SIN 

— 

DSORTF 

EL? 

els 
EL 1 5 
EL 1 6 

SORT 

— 
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VI. BRIEF DESCRIPTION OF SUBPROGRAMS 


An alphabetical and numerical index, cross-referencing the 124 subprograms 
contained in RESIST, appears in Figure 6-1. A brief description of each sub- 
program is also presented in the order in which it appears in the FORTRAN 
program. 
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Subprogram Order* Subprogram Order* Subprogram Order* Subprogram Order* Subprogram Order* 


ALPHABETICAL LISTING OF SUBROUTINES 


ADD 

BIGBRD 

COMPAR 

DFIND 

DSQRTF 

ELI 5 

ENTMPY 

FTEMP 

GETROW 

^INPUT 

KPR I NT 

LOCDEF 

MEMBIN 

MSHWRT 

NSTAK 

PODSYM 

PROOFS 

OBAC 

OFACT 

RDM 

SAMAIN 

SLOAD 

SREVN2 

TTNTER 

UNPACK 



ALARM 

05 6 

ANDOR 

062 

ARACE 

050 

ASTACK 

BOUND 

00 8 

BPMASS 

055 

BRICK 

006 

CARDIN 

COPY 

02* 

DABSF 

121 

DCLOSE 

11? 

DCOSF 

DLGRN 

022 

DOTPRO 

061 

DREAD 

107 

DSTNF 

DSRAN 

OB 4 

DWRITE 

104 

EDGMAS 

049 

FLAS 

EL 1 6 

04 6 

EL 2 

04 3 

ELS 

044 

ENMMPY 

EOFOl 

117 

FLXTSL 

048 

FMSTS 

026 

FSAKG 

GAUSS 

02 8 

GEOBC 

005 

GET 

108 

GETDIM 

HEXEL 

02 9 

HE XT IL 

082 

HOTDOT 

0_9_8_ 

ICMULT — 

— I-N-T-EMP — 
LAGR AN 

-0-1-8- 

01 7 

— ISTAK 

LDB 

030 

115 

JACOB 

LGRANG 

031 

023 

KMOSO 

LOADIN 

LOWTRT 

069 

MAIN 

001 

MATB 

09? 

MATL 

MEMGEN 

01 0 

MESH 

012 

MPR INT 

097 

MPROP 

MSOUT 

095 

MSTACK 

054 

MSTAK 

032 

MULT 

ORTHOG 

080 

PACK 

113 

PARTIN 

034 

PLB 

PREAD 

106 

PREP 

059 

PRINT3 

051 

PRITE 

PSTRES 

094 

PUT 

105 

PUTLAB 

099 

PUTROW 

OB SOL 

076 

OCHOL 

071 

ODLD 

0 35 

ODST 

OFOR 

075 

OF SOL 

074 

OPASS 

073 

OSVEC2 

READMT 

067 

REND 

110 

RESTOR 

064 

REVERS 

SB MAIN 

052 

SCAPRO 

057 

SETUP 

003 

SHAPE 

SPLITS 

04 2 

SR AIN 

085 

SRAIN2 

086 

SRA1N3 

STURM 

05 8 

SUB 

091 

TEMPER 

020 

TIME A 

TITLES 

118 

TLGRN 

021 

TRAN 

079 

UNDCUT 

WEND 

109 

WRAKG 

039 

WRDSP 

088 



MAIN 

00 1 

SAMAIN 

BRICK 

00 6 

UNDCUT 

LOADIN 

01 1 

MESH 

T INTER 

016 

LAGRAN 

TLGRN 

021 

DLGRN 

FMSTS 

026 

FSAKG 

JACOB 

031 

MSTAK 

QDST 

036 

SHAPE 

INPUT 

04 1 

SPLITS 

ELI 6 

04 6 

SREVN2 

PR I NT3 

05 1 

SBMAIN 

ALARM 

056 

SCAPRO 

DOTPRO 

06 1 

ANDOR 

ENMMPY 

066 

READMT 

OCHOL 

071 

KPR I NT 

QfeSOL 

076 

QB AC 

LOCDEF 

081 

HEXTIL 

SRAIN2 

08 6 

SRATN3 

SUB 

091 

MATB 

COMPAR 

096 

MPR I NT 

GETDIM 

10 1 

GETROW 

PREAD 

106 

DREAD 

DFIND 

11 1 

DCLOSE 

PLB 

1 1 6 

EOFOl 

DABSF 

12 1 

DCOSF 



SETUP 

003 

CARDIN 

BOUND 

008 

MEMBIN 

MSHWRT 

013 

MPROP 

INTEMP j 

018 

FTEMP 

LGRANG 

023 

COPY 

GAUSS 

028 

HEXEL 

NSTAK 

033 

PARTIN 

SLOAD 

038 

WRAKG 

EL 2 

043 

EL5 

FLXTBL 

048 

FDGMAS 

ASTACK 

053 

MSTACK 

STURM 

058 

PREP 

RDM 

063 

RESTOR 

ICMULT 

068 

LOWTRI 

OPASS 

073 

OF SOL 

REVFRS 

078 

TRAN 

KMOSO 

083 

DSRAN 

WRDSP 

088 

PODSYM 

FNTMPY 

093 

PSTRES 

HOTDOT 

098 

PUTLAB 

PRITE 

103 

DWRITE 

GET 

108 

WEND 

PACK 

113 

UNPACK 

TITLES 

118 

BIGBRD 

DSTNF 

123 

DSQRTF 



GFOBC 

*EWGEN 

MATL 

TEMPER 

ETAS 

ISTAK 

COLD 

PROCES 

FLI5 

APACE 

BPMASS 

GSVEC2 

HJLT 

OF ACT 

OF OR 

CRTHOG 

SPAIN 

ADD 

MSOUT 

PUTROW 

PUT 

REND 

LDB 

TI ME A 





*Order in which subprogram appears in program 


Fig. 6-1 Subprogram Nomenclature and Numerical Index 
































' DESCRIPTION OF SUBROUTINES 


1 MAIN - Initializes input/output variables. 

2 SAMAIN - Supervisor that calls all input data processing and iteration 
routines. 

3 SETUP - Assigns input/output unit and file numbers. 

4 CARDIN - Prints problem options and calls data generation and pre-processing 
routines. . 

5 GEOBC - Writes primary structure coordinates and boundary condition on data 
set. 

6 BRICK - Generates primary structure coordinates, corresponding to brick tile 
configuration. 

7 UNDCUT - Generates primary structure coordinates corresponding to undercut 
tile configuration. 

8 BOUND - Generates primary structure boundary conditions. 

9 MEMBIN - Picks out nodes corresponding to each primary structure element. 

10 MEMSEN - Stores all material property information and nodes for each primary 
structure element on the data set. 

11 LOADIN - Reads in primary structure mechanical and thermal loading data (for 
statics problems only). Generates applied load matrix. 

12 MESH - Generates nodal geometry for brick and undercut tiles. Prints co- 
ordinates, element and node maps if requested. 

13 MSHWRT - Generates material properties for the elements and stores them, 
along with geometry, on the data set. 

14 MPROP - Reads in material property data for primary structure elements, tile 
arrestor, isolator and RSI and prints them. 

15 MATL - Stores material property data for tiles in matrix form to be used in 
stiffness matrix calculations. 
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16 TINTER - Calls the Lagrangian interpolation routine for temperature 
dependent material properties of tiles. 

17 LAGRAN - Performs Lagrangrian interpolation of tile material property tables. 

18 INTEMP - Reads in data for the static thermal loading of tiles via one of 
three options - uniform temperature, Lagrangian interpolation or element node 
temperature . 

19 FTEMP - Determines tile nodal temperatures. 

20 TEMPER - Determines whether or not x, y, z coordinate arguments are located 
wi-thin the -interpolation~tables". 

21 TLGRN - Performs three dimensional Lagrangian interpolation. 

22 DLGRN - Performs two dimensional Lagrangian interpolation. 

23 LGRANG - Performs single argument Lagrangian interpolation. 

24 COPY - Copies a matrix from one data set to another. 

25 ELAS - Supervises the tile stiffness matrix generation. 

26 FMSTS - Forms the strain matrix for the tile element centroid. 

27 FSAKG - Supervises the formation and stacking of the tile element stiffness 

matrices . 

28 GAUSS - Forms Gaussian integration table for 2 thru 6 points. 

29 HEXEL - Generates tile element stiffness and load matrices and supervises 

their stacking. 

30 ISTAK - Writes tile element stiffness and load matrices with their stacking 
indices on data set. 

31 JACOB - Forms Jacobian and its inverse. 

32 MSTAK - Returns stacking index for tile element stiffness matrix. 

33 NSTAK - Stacks total tile stiffness matrix and load vector. 

34 PARTIN - Partitions out specified rows of tile stiffness or load matrices. 

35 QDLD - Forms quadrilateral (coating) load matrix. 

36 QPST - Forms quadrilateral (coating) element stiffness matrix. 

37 SHAPE - Contains the hexahedron shape function. 
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38 SLOAD - Writes tile load vector if requested. 

39 WRAKG - Prints tile element stiffness matrices if requested. 

40 PROCES - Keeps a running total on the number of primary structure degrees of 
freedom at each node and prints them out, along with nodal coordinates. 

41 INPUT - Calls the appropriate finite element subroutine to calculate individual 
primary structure element stiffness matrices. Also used for mass matrices in 
vibration problems. 

42 SPLITS - Computes stacking indices and band-widths for primary structure 
element stiffness and mass matrices. 

43 EL2 - Generates beam finite element stiffness and mass matrices in global 
coordinates . 

44 EL 5 - Generates global stiffness and TPS mass matrices for a planar isotropic 
quadrilateral membrane finite element. 

45 EL15 - Generates stiffness matrix in global coordinates for triangular plate 
bending finite element. 

46 ELl 6 - Generates element stiffness and mass matrices for qualrilateral 
bending element. It calls ELI 5 four times to assemble four triangles into a 
"quad" with interior degrees of freedom included which are subsequently con- 
densed out of local stiffness matrix before transformation to global co- 
ordinates. 

47 SREVN2 - Obtains inverse of small order matrix through Gauss-Jordan elimina- 
tion scheme with partial pivoting and stores it in the position of the 
original matrix. 

48 FLXIBL - When the primary structure has flexible boundary conditions, this 
routine identifies the appropriate boundary elements and modifies the cor- 
responding finite element stiffness-matrices to include the effects of the 
flexible boundaries. 

49 EDGMAS - Adds overhung rotatory mass item into beam mass matrix of elements 
on A and B edge of primary structure. 

50 ARACE - Eliminates rows and/or columns of primary structure element stiff- 
ness or mass matrices, according to boundary conditions associated with 
element. 
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51 PRINT3 - Prints all primary structure stiffness and mass matrices, - if 
requested. 

52 SBMAIN - Controls whole iterative procedure of the program including: 
primary structure matrix stacking, equation solution routines, approximate 
frequency and mode shape calculations, convergence criteria tests and final 
printout of stresses. 

53 ASTACK - Reads element stiffness matrices and associated stacking indices 
from data set and stacks them into total primary structure stiffness matrix. 

54 MSTACK - Reads element mas s matrices and associated stackin g indices_from 

data set and stacks them into total primary structure mass matrix. 

55 BPMASS - Reads element mass matrices and associated stacking indices from 
data set and stacks only those associated with plates and beams into total 
primary structure diagonal mass matrix. 

56 ALARM - Automatically reduces size of a vibration problem and computes 
approximate frequencies and mode shapes. 

57 SCAPRO - Accumulates products of row and column elements of matrices. 

58 STURM - Applies bisection technique to obtain eignvalues of a real symmetric 
tri-diagonal matrix. 

59 PREP - Provides routine STURM with critical information for location and 
number of roots. 

60 QSVEC2 - Computes unit eigenvector of matrix, given diagonal and off-diagonal 
entries of this Householder tri-diagonal matrix and a good approximate root. 

61 DOTPRO - Calculates "dot" product of - two vectors. 

62 ANLOR - Logical "and" function and logical "or" function. 

63 RDM - Generates random numbers uniformly distributed between 0 and 1. 

64 RESTOR - Generates symmetric matrix from its lower triangle. 

65 MULT - Controls manipulation of two matrices on data sets which are to be 
multiplied together in packed form and stores the result on a third data set. 

66 ENMMPY - Multiplies two matrices stored in core to produce partial product of 
two larger matrices which do' not fit in core. 
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67 READMT - Reads column vector from data set into core. 

68 ICMULT - Multiplies adjacent elements of two column vectors and 
sums results to form scalar quantity. 

69 LOWTRI - Generates lower triangle from full matrix stored on data set. 

70 QFACT - Reads positive definite symmetric stiffness matrix from data set and 
sets up indices for subroutine QCHOL to get lower triangle of an LL^ 
decomposition one row at a time. 

71 QCHOL - Performs Cholesky factorization of stiffness matrix. 

72 KFRINT - Prints non-zero elements of lower triangular primary structure 
stiffness matrix. 

73 QPASS - Reads over required number of rows of lower triangle decomposition 
of stiffness matrix to aid in zoning for forward solution of system of 
equations; 

74 QFSOL - Sets up and manages zoning for forward solution of L-^ X = L"-4f = Z, 
using L and Y and getting Z. 

75 QFOR - Computes forward solution of L - ^ Y = Z. 

76 QBSOL - Sets up and manages zoning for backward solution of L 1 X = Z to 
obtain X. 

77 QBAC - Computes backward solution of L T X = Z. 

78 REVERS - Reverses order of rows of matrix on data set when matrix will not 
fit in core. 

79 TRAN - Computes transpose of matrix which is stored on data set and will not 
fit in core. 

80 ORTHOG - Normalizes elements of column vector. 

81 LOCDEF - Determines primary structure nodes and deflections associated with 
boundary nodes of each tile and generates new global load vector upon comple- 
tion of tile analysis. 

82 HEXTIL - Controls partitioning of tile stiffness and load matrices, solution of 
equations, computation of tile portion of Rayleigh Quotient and calculation 

of tile stresses and strains. 
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83 KMOSQ, - Computes K - u> M for tiles and primary structure. 

84 DSRAN - Calculates tile member stresses and strains. 

85 SRAIN - Computes and prints stresses and direct strains for each tile. 

86 SRA.IN 2 - Computes and prints stresses for each tile isolator and arrestor. 

87 SRAM 3 - Computes- and prints stresses and total strains for coating for 
each tile. 

88 WRDSP - Prints TPS displacements for each tile. 

89 PODSYM - Calls equation sol ving_r.outines-in-appronr-iate-order-i 

90 ADD - Adds two matrices stored on data sets which may not fit in core and 
stores result on third data set. 

91 SUB - Subtracts two matrices stored on data sets which may not fit in core 
and stores result on third data set. 

92 MATB - Controls manipulation of matrices A and B which are read from data 
sets, zones arrays to calculate A X B and stores result on third data set. 

93 ENTMPY - Multiplies two matrices stored in core to produce partial 
product A^B of two larger matrices A and B which do not fit in core. 

94 PSTRES - Calculates and prints plate strains and stresses if requested. 

95 MSOUT - Prints primary structure deflections at each node. 

96 COMPAR - Computes maximum deflection, minimum deflection difference and 
maximum convergence parameter between iterations. 

97 MPRIET - Prints non-zero elements of lower-triangular primary structure mass 
matrix. 

98 HOTDOT - Forms inner product of two vectors. 

99 PUTLAB - Puts label on data set containing matrix dimensions and prints out 
label information, if requested. 

100 FUTROW - Writes row of matrix on data set in format designated by packing 
factor. 

101 GETDIM - Gets label of matrix from data set and prints out label information, 
if requested. 
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102 GETROW - Reads row of matrix from data set in format designated by packing 
factor. 

103 ERITE - Converts number of bytes in array into words by dividing by four. 

10b DWRITE - Converts number of bytes in array into words by dividing by four. 

105 PUT - Writes singly dimensioned array onto data set. 

106 FREAD - Converts number of bytes in array into words by dividing by four. 

107 DREAD - Converts number of bytes in array into words by dividing by four. 

108 GET - Reads singly dimensioned array from data set. 

109 WEND - Called at end of a series of PRITE'S. 

110 REND - Called at end of a series of PREAD'S. 

111 DFIND - Rewinds data set to appropriate file. 

112 DCLOSE - Rewinds data set. 

113 PACK - Packs rows of matrix so they may be written on a data set in an effi 
cient manner. This is done by representing strings of zeroes by a single 
fixed point negative integer where the value of the integer represents the 
number of zeroes in the string. Non-zero numbers are preceded by a fixed 
point number indicating the number of non-zero numbers that follow. 

114 UNPACK - Unpacks rows of matrices that have been packed by subroutine PACK. 

115 LDB - Reads and lists input data from input stream and generates data file 

for user's program. 

116 PLB - Starts a new page and/or skips a number of lines. 

117 E0F01 - Checks for an end of file on input data set unit, sets end of file 

control word option to its proper value and reads a two card problem title. 

118 TITLES - Prints title at top of page. 

119 BIGBRD - Prints RESIST title sheet before program output. 

120 TIMEA - Calculates and lists computer time at desired intervals within 
program. 

121 DABSF - Takes absolute value of floating point number. 
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122 DCOSF - Computes cosine of floating point argument. 

123 DSINF - Computes sine of floating point argument. 

124 DSQRTF - Computes square root of positive' floating point number. 
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VII. TYPICAL SUBPROGRAM OUTPUT 

The following listings display subprogram output for two problem types 
with all the debugging clues turned on. One execution is for a statics option 
problem and the second is for a dynamics option problem. The name of the parti- 
cular subprogram which generated each set of output is printed in bold type on 
the rotated right hand side of each photographed computer page. 
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THE ICWFR TRIANGLE OF THE TOTAL STIFFNESS MATRIX OF THE PRIMARY STRUCTURE 
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MATRIX (MASS S ) 9 ROWS* 9 COLUMNS, FROM PROGRAM EL2 
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the lower triangle of the total stiffness matrix of THE primary structure 
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EL MASS * I/O UNIT = 8« FILE = 1. ROWS = 8, COLUMNS = 24 



THE LOWER TRIANGLE OF THE TOTAL MASS MATRIX OF THE PRIMARY STRUCTURE 
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DESIRED MODES = 3 

REORTHOGONALI ZATI ONS 
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THE LOWER TRIANGLE OF THE TOTAL STIFFNESS MATRIX OF THE PRIMARY STRUCTURE 
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